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We present a wide study on the comparison of different shadowing models and their influence 
on J/i/> production. We have taken into account the possibility of different partonic processes 
f~| . for the cc-pair production. We notice that the effect of shadowing corrections on J/tp pro- 

^ l' duction clearly depends on the partonic process considered. Our results are compared to the 

available data on dAu collisions at RfflC energies. We try different break up cross section for 
each of the studied shadowing models. 



1 Introduction 

> 

' The interest devoted to J/i/j production has not decreased for the last thirty years. It is motivated 

Q^ ■ by the search of the transition from hadronic matter to a deconfined state of QCD matter, the 

i so-called Quark-GIuon Plasma (QGP). The high density of gluons in the QGP is expected to 

I hinder the formation of quarkonium systems, by a process analogous to Debye screening of the 

' electromagnetic field in a plasma^]. 

■ The results on J ftp production measured by the PHENIX experiment at the BNL Relativistic 

Heavy Ion Collider (RHIC) show a significant suppression of the J/V' yield in AuAu collisions 
^ I at y^ijv7V=200 GeV . Nevertheless, PHENIX data on dAu collisions!^ have also pointed out 

^ ■ that Cold Nuclear Matter (CNM) effects play an essential role at these energies. 



All this reveals that, in fact, the interpretation of the results obtained in nucleus- nucleus 
collisions relies on a good understanding and a proper subtraction of the CNM effects, known 
to impact the J/ip production in proton (deuteron)-nucleus collisions, where the deconfinement 
can not be reached. 

It is our purpose here to develop an exhaustive study of these effects. This includes shado- 
wing, i.e. the modification of the parton distribution of a nucleon in a nucleus, and final-state 
nuclear absorption. 

Moreover, we have recently noticed^ that the impact of gluon shadowing on J/tp production 
does depend on the partonic process producing the cc and then the J/ip. Because of this, we have 
included here two different production mechanism: g + g^cc^J/ij) that corresponds 

to the picture of the Colour Evaporation Model (CEM) at LO, and g + g ^ J/tp -\- g, as in the 
Color Singlet Model (CSM), but including the s-channel cut contributions. 

In practice, we have proceeded as follows: we have developed a Glauber Monte-Carlo code 
where we have interfaced the two different partonic processes for cc with the CNM effects, 
namely the shadowing and nuclear absorption, in order to get the J/ip production cross sections 
for pA and AA collisions. We have considered four different models for the shadowing effects. 
We shall compare our results with the experimental measurements on dAu collisions presently 
available at RHIC. 



2 The Model 



To describe the J /ip production in nucleus collisions, our Monte-Carlo framework is based on the 
probabilistic Glauber model, the nuclear density profiles being defined with the Woods-Saxon 
parameterisation for any nucleus A > 2 and the Hulthen wavefunction for the deuteron. The 
nucleon-nucleon inelastic cross section at ^snn = 200 GeV is taken to ajyN = 42 mb and the 
average nucleon density to po = 0.17 nucleons/fm'^. For each event (for each AB collision) at a 
random impact parameter b, the Glauber Monte-Carlo model allows us to determine the number 
of nucleons in the path of each incoming nucleon, therefore allowing us to easily derive the 
number A'coU of nucleon-nucleon collisions and the total number A^part of nucleons participating 
into the collision. In order to study the J/tp production, we need to implement in our Monte 
Carlo the following ingredients: the partonic process for the cc production and the CNM effects. 



2.1 Partonic process for the cc production 

Most of the studies of J/^jJ production rely on the assumption that the cc pair is produced by 
the fusion of two gluons carrying some intrinsic transverse momentum kx- The partonic process 
being a 2 ^ 1 scattering, the sum of the gluon intrinsic transverse momentum is transferred to 
the cc pair, thus to the J/ip since the soft hadronisation process does not modify the kinemat- 
ics. This corresponds to the picture of the Colour Evaporation Model (CEM) at LO. In such 
approaches, the transverse momentum of the J/tp entirely comes from the intrinsic transverse 
momentum of the initial gluons. 

However, such an effect is not sufficient to describe the Pt spectrum of quarkonia produced 
in hadron collisions^. Most of the transverse momentum should have an extrinsic origin, i.e. the 
J/ip^s Pt would be balanced by the emission of a recoiling particle in the final state. The J ftp 
would then be produced by gluon fusion in a 2 ^ 2 process with emission of a hard final-state 
gluon. This emission, which is anyhow mandatory to conserve C-parity, has a definite influence 
on the kinematics of the J/ip production. Indeed, for a given J/ip momentum (thus for fixed y 
and Pt), the processes discussed above, i.e. g + g ^ cc ^ J/tp {+X) and g + g ^ J /ip + g, will 
proceed on the average from gluons with different Bjorken-x. Therefore, they will be affected 
by different shadowing corrections. From now on, we will refer to the former scenario as the 
intrinsic scheme, and to the latter as the extrinsic scheme. 

In the intrinsic scheme, the initial gluons carry a non-zero intrinsic transverse momentum, 
which is transferred to the J/ip . Note that, following , we do not neglect the value of the J/\p^s 
Pt in this simplified kinematics. In fact, in this scheme, we use the fits to the y and Pt spectra 
measured by PHENIX in pp collisions at y'ijvTV = 200 GeV as inputs of the Monte-Carlo. 
The measurement of the J/ ip momentum completely fixes the longitudinal momentum fraction 
carried by the initial partons: 

xi,2 = exp (±y) = x%{y, Pt), (1) 



with the transverse mass my = y M'^ + Pj,, M being the J/V' mass. 

On the other hand, in the extrinsic scheme, information from the data alone - the y and Pt 
spectra - is not sufficient to determine xi and X2. Indeed, the presence of a final-state gluon 
authorizes much more freedom to choose {xi,X2) for a given set {y,PT). The four-momentum 
conservation explicitely results in a more complex expression of X2 as a function of {xi,y,PT): 

_ ximTy/SNNe~y - 
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Equivalently, a similar expression can be written for xi as a function of {x2,y, Pt)- Even if the 
kinematics determines the physical phase space, models are anyhow mandatory to compute the 



proper weighting of each kinematically allowed (xi,X2). This weight is simply the differential 
cross section at the partonic level times the gluon Parton Distribution Functions (PDFs), i.e. 
g{xi, fif)g{x2, ^J-f) dagg^j/^^g/dy dPT dxidx2. In the present implementation of our code, we 
are able to use the partonic differential cross section computed from any theoretical approach. 
For now, we use the one from^ which takes into account the s-channel cut contributions^" to the 
basic Color Singlet Model (CSM) and satisfactorily describes the data down to very low Pt . 



2. 2 Shadowing 

To get the J /ip yield in pA and AA collisions, a shadowing-correction factor has to be applied to 
the J /ip yield obtained from the simple superposition of the equivalent number of pp collisions. 
This shadowing factor can be expressed in terms of the ratios Rf" of the nuclear Parton Distri- 
bution Functions (nPDF) in a nucleon of a nucleus A to the PDF in the free nucleon. We will 
consider three different shadowing models for comparation: EKS98 and nDSg^ 

at LO. 

These models provide the nuclear ratios at a given initial value of Qq which is assumed 
large enough for perturbative evolution to be applied: Qg = 2.25 GeV^ for EKS98, = 
1.69 GeV^ for EPS08 and Ql = 0.4 GeV2(0.26 GeV^) for nDSg and perform their evolution 
through the DGLAP evolution equations to LO accuracy in the case of EKS98 and EPS08 and 
to LO and NLO accuracy in the case of nDSg. The spatial dependence of the shadowing is not 
given in the above models. However, it has been included in our approach, assuming that the 
inhomogeneous shadowing is proportional to the local density^^. 

The nuclear ratios of the PDFs are then expressed by: 



The numerical parameterisation of Rf{x,Q ) is given for all parton flavours. Here, we restrain 
our study to gluons since, at high energy, J/V' is essentially produced through gluon fusion^. 



2.3 The nuclear absorption 

The second CNM effect that we are going to take into account concerns the nuclear absorption. 
In the framework of the probabilistic Glauber model, this effect refers to the probability for the 
pre-resonant cc pair to survive to the propagation through the nuclear medium and is usually 
parametrised by introducing an effective absorption cross section (Tabs- It is our purpose here 
to compare different absorptive a within the two partonic cc production mechanisms -intrinsic 
and extrinsic- and for the three shadowing models cited above. 



3 Results 



In the following, we present our results for the J/^l^ nuclear modification factor: 

dN'ji^Q {dNpp^) is the J/'^ yield observed in AB {pp) collisions and (Ncou) is the average number 
of nucleon-nucleon collisions occurring in one AB collision. In the absence of nuclear effects, 
Rab should equal unity. 

We will restrict ourselves to dAu collisions, since only CNM matter effects are at play here, 
so they provide the best field for the study of the shadowing and the nuclear absorption. 



We have used PHENIX measurements of -R^Au in order to compare the different shadowing 
models. In Fig. [H we have computed our results in the different shadowing frameworks for four 
(Tabs for each of the shadowing models considered in both the intrinsic and extrinsic scheme. 

We have also evalute the best fit for dabs, following the method used by PHENIX in"^ and 
By using the data on RdAu versus rapidity, we have obtained the best foi^ the EPS08 
model, computed in the extrinsic scheme for a (Tabs = 3.6 mb. 




Figure 1: J/tp nuclear moiiification factor in dAu collisions at y^sjvjv = 200 GeV versus y for the differ- 
ent sliadowing models in the intrinsic (up) and extrinsic (down) scheme. For complementary figures see: 
http : //phenix-f ranee . in2p3.fr/ sof tware/ jin/ index . html ^ 
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